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ABSTRACT 

Electrons in a spherical ultracold quasineutral plasma at temperature in the Kelvin 
range can be created by laser excitation of an ultra-cold laser cooled atomic cloud. 
The dynamical behavior of the electrons is similar to the one described by conven- 
tional models of stars clusters dynamics. The single mass component, the spherical 
symmetry and no stars evolution are here accurate assumptions. The analog of binary 
stars formations in the cluster case is three-body recombination in Rydberg atoms 
in the plasma case with the same Heggie's law: soft binaries get softer and hard bi- 
naries get harder. We demonstrate that the evolution of such an ultracold plasma is 
dominated by Fokker-Planck kinetics equations formally identical to the ones control- 
ling the evolution of a stars cluster. The Virial theorem leads to a link between the 
plasma temperature and the ions and electrons numbers. The Fokker-Planck equation 
is approximate using gaseous and fluid models. We found that the electrons are in a 
Kramers- Michie-King's type quasi-equilibrium distribution as stars in clusters. Know- 
ing the electron distribution and using forced fast electron extraction we are able to 
determine the plasma temperature knowing the trapping potential depth. 

Key words: stellar dynamics - plasmas - atomic processes - (stars:) binaries: general 



1 INTRODUCTION 

One challenge of astrophysics is to unde rstand the dynamics 
of glob ular star clusters (for a review see lMevlan and Heggid 
il997l) ') because they are test systems for dynamical theoret- 
ical models such as N-Body modeling, Monte Carlo simula- 
tio ns of Fokker-Planck equation s , gas model, scaling m odels, 
... jBinnev and TremainjllQSTt iLvman Spitzejll987ft . The 
evolution of globular clusters is dominated by two (or three)- 
body relaxation, evaporation of stars, tidal truncation and 
stars evolution. This leads to a very complex evolution and 
the theoretical models are thus simplified with many approx- 
imations sometimes far from reality. Furthermore it is not 
possible to observe the evolution of a given cluster because 
observation gives only an instantaneous picture. 

In this letter we propose to study a real system which 
corresponds to the most usual assumptions used in conven- 
tional models of stars clusters: a single mass component, an 
almost perfect spherical symmetry and no stars evolution. 
This system which can then be efficiently compared with 
theory, namely, an ultra-cold plasma is realized, controlled 
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and studied in a small laboratory. An ultra-cold plasma can 
be formed by laser excitation of an ultra-cold (T ~ 100 nK) 
atom ic sample and has been first realized bv lKillian et alJ 
I1999I The physics of ultracold plasmas have strong similari- 
ties with the physics of globular stars clusters. Both systems 
are spherically symmetric, radially limited due to tidal forces 
for clusters or due to applications of an external magnetic or 
electric field for plasma. The key point, developed through 
the whole article, is that both systems are driven by the same 
kinetic equations. Indeed, they are subject to the same in- 
verse square type forces if one uses, for the plasma case, a 
new strong negative "gravitational constant" G' defined by: 



G' 
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-2.78 X lO^V^kg" 



(1) 



Analog of binary stars are excited Rydberg atoms and three 
body recombinations play the same role in both systems. 
For instance in a cluster the energy source, in post-collapse 
evolution, is provided by binaries formation. Similarly in 
expanding plasma an heating is provided by three body re- 
combination and Rydberg atoms formation. However, some 
aspects are not perfectly matched in both systems. For in- 
stance a globular cluster orbits around its host galaxy and 
is therefore submitted to centrifugal and tidal forces. On 
the contrary, the plasma is not orbiting. Nevertheless, ex- 
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ternal field (electric or magnetic) might mimic such forces. 
In fact the main goal of this article is to demonstrate that, 
in ultra-cold plasma, the electrons dynamics confined by the 
ionic potential is almost identical to equal mass non evolving 
stars dynamics in a globular cluster. 

The analogy has already been u sed in A''-body simu- 
lation by iKuzmin and O'Neill l2002bllal who simulated the 
plasma behavior using a modified Aarseth's code usually 
devoted to cluster studies. In fact the modifications were 
severe (S. V. Kuzmin personal communication) and it is 
not possible to simply use globular clusters N-Body system 
equations to fully study ultra-cold plasma system using only 
physical constant replacement. It is nevertheless possible to 
precisely study electrons (mass rrie, charge Qe < 0, spatial 
density rie) behavior within the ionic external potential (ion 
mass rrii, ion charge Qi = —qe > 0, spatial density rii) using 
a tr ue analogy with a two mass component globular clus- 
ter. IVanhaecke et al]l2004l were the first to notice a formal 
analogy and have developed a close analogy using a low- 
ered Maxwellian King's type distribution for the electrons 
in the plasma. Similar distribution has then been used by 
FPohl et al. (.2004a.bi) . This letter is devoted to an extension 
of the analogy. The article has been written to be accessible 
to the ultra-cold plasma community as well as the stellar 
dynamics community. Lots of ideas, equations or models are 
not fully solved in this paper. Nevertheless we have though 
that it is fair and interesting to present preliminary equa- 
tions in order to stimulate further works. The outline of the 
paper is as follows. In section |21 we discuss the analogy be- 
tween the plasma and the cluster system through scaling 
relation in mass, length and times units. In section |3| we 
present the basics of ultra-cold plasma physics and the role 
of the three body mechanism. In section 0] we show that, 
the Fokker-Planck equation is exactly identical for plasmas 
and clusters system. In section|S]we briefly study the evolu- 
tion of the quasi-equilibrium distribution using gaseous and 
fluid model for ions and electrons. Section |S| is devoted to 
the orbit average Fokker-Planck equation. We discuss some 
physical consequence as electron evaporation rate. We detail 
also some approximate solution as a King's type equilibrium 
distribution for electrons in the plasma in analogy with glob- 
ular star cluster dynamics. In section |7| we detail some ex- 
perimental results. For instance we experimentally extract 
electrons trapped in an ionic cloud by an external electric 
fleld to estimate the electron temperature of the plasma. We 
finally summarize our results and their implications for both 
astrophysics and plasma physics. 



2 DYNAMICAL ANALOGY THROUGH 
SCALE UNITS 

In table^we give the typical parameters of both systems. It 
is then obvious that the natural units (mass, size, crossing 
time) for both system are severals tens orders of magnitudes 
different. But, if we use the units nip = ms,rp = 200 fim, 
tp = 100ns for plasma and Mc = Mq,Rc = 1 pc, Tc = 
10 Myrs for clusters, the systems looks quite similar. Indeed 
the fc"^ electron evolution in the plasma is described by the 
Newton's law, with dimensionless scaled notations f — r/rp, 



t = t/tp and rUe 
me — — ~ 

dt^ 



■ me /nip: 
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where Upfii = Ui is the ion density. 

By comparison the full dynamics for an Mi mass star 
in a two component cluster (stars mass Afi and M2 only) 
in uppercase dimensionless units R — R/ Rc, T — t/Tc and 
Ml = Ml /Mc is given by: 
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providing f — R, rh — AI and t = T the equation are ex- 
actly identical with similar coupling constant; G'm^ — ^ 



and GMl -j^f^ on the order of unity (see equation JTJ and 
values in table 0. 

However it is not possible to simply use globular clus- 
ters N-Body code to study ultra-cold plasma system. This 
is due to the poorly behavior of the second type of particles; 
ions or Kh- A first solution is to mod ify an existing code 
devote d to cluster studies as done by iKuzmin and O'Neill 
i2Qfl2^0. In a similar manner IPohl et al.ll2004bl have used 
a molecular dynamics simulation with ionic correlations 
based on the treecode originally designed for astrophysics 
problems by J. E. Barnes. Another strategy is to use a 
Monte-Carlo method. Detail of one Monte Carlo method 
used for ultra-cold plasma with Three-Body Recombination 
(TBR) process and radiative atomic lifetime is given by 
iRobicheaux and Hansonl2003l . Analogy with cluster code as 
discussed in the MODEST group (Modeling DEnse STel- 
lar systems, http://www.manyb0dy.0rg/m0dest/l is obvi- 
ous and this relation would be probably fruitful for further 
studies. 



3 ULTRACOLD PLASMA 
3.1 Ultracold plasma formation 

The ultracold neutral plasma physics experimentally be- 
gan at NIST in 1999 with laser photoionization of a laser- 
cooled (in the microkelvin or millikelvin range) sample. 
Since, several groups have done similar work with trapped 



met astable xenon, cesium or rubidium atomsffor a review 
see JGallaeher et all I2OO3I : iKillian et alJ 12003)1. Recently 
strontium atoms have been also used bv lsimien et alJl2004l . 
This is a promising experiment because the strontium ion 
can be further laser cooled. This might be a way to controlled 
the ion motions. A typical experimental setup is presented 
in figure^ Due to the small electron over ion mass ratio, the 
electrons have an initial kinetic energy kj = SknTe /2 (ks is 
the Boltzmann constant) which is, at first glance, believed to 
be approximately equal to the difference between the photon 
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Table 1. Typical units value for globular system and for ultra-cold plasma. 





Stars cluster 


Electrons in Plasma 


Number of particules 


5 X 10^ 


10'' 


System Mass 


5 X 10'' Mq lO'^'" kg 


10" me Si 9 X 10"^^ kg 


Size (Core radius) 


1 pc « 3 X 10^'' in 


200 /im = 2 X 10"" m 


Peak density 


8 X lO^M0pc-^ 


10"' cm-^ 


Crossing time 


10 Myrs Ri 3 X lO^'' s 


100 ns = 10"''' s 


Velocity dispersion a^, 


7 km/s 


50 km/s 


Relaxation time 


10 Gyrs 


10 ns 


Coulomb logarithm : In A 


10 


4 



energy and the ionization threshold. T2 is the photoionisa- 
tion electronic temperature. Parameters can easily be exper- 
imentally tuned in the following range: < 10* electrons 
at temperature T2 < 1000 K embedded in Ni < 10* ions 
(initial temperature ~ 1 mK). As indicated in figure Q 
the spherical gaussian atomic sample is ionized by a gaussian 
laser with an intensity profile given by / = Ioc~ i-z ; /m 
This leads to a cylindrical initial plasma shape with a gaus- 
sian radius along the longitudinal [x) laser propagation axe, 

a[t — Q) = ao ~ 250 /xm, and with yj^^^ gaussian radius 

along the radial [y, z) axes. Experiments are sometimes done 
with w « (T so the spherical symmetry is not at all perfect. In 
the following we will assume ti; ^ cr to restore the gaussian 
spherical symmetry. We then have ni{r, t) = n1e~^ ^'^^'^ 
with density n° < 10^^ cm~"^. With typical initial value 
Te = 50 K, ni ^ rii — l(f cm~'^, the plasma parameters 
are the following: ion (or electrons) number Ni = Ne ~ 
250 000, Debye screening length Ao = eoksTe / (qinl) ^ 
15 ^m, the Wigner-Seitz radius (interparticule spacing) 
dws = (47rn°/3)~^'''' ~ 6 /im. Landau length tl ~ 
/(47reofcs7e) « 0.3 /^m and the Coulomb logarithm 
InA = ln(2AB/ri,) = \n{^-Kq:''{eokBTef^^{nl)-^'^) ^ 
4. There is several definitions for the electron relaxation 
time with slightly differen t numerical factor (as for cluster 
iLouis and Spurzemlll99lh ). We choose here the relaxation 
times te to be defined by te = 9i'J/(16x/7rG'^mene In A) ^ 
15 ns for electron with velocity Ve, chosen for the numerical 
results at its r.m.s. value = {{ve — {ve)Y)^^'^ ~ 50 km/s. 
The relaxation time is fast enough to have electrons al- 
ways in quasi-equilibrium. On the contrary the ions are al- 
most never thermalized with the electrons during the typ- 
ical 100 /xs plasma expansion lifetime. In the previous def- 
initions the temperature has to be understood as the radi- 
ally dependent thermodynamic (velocity average) temper- 
ature imecrj; = ^ksTeir). It is worth to note that the 
plasma is a kinetic plasma because the thermal energy is 
higher than the Coulomb interaction energy {rL < aws)- 
In fact one major goal of the ultra-cold plasma commu- 
nity is to reach the opposite situation, the so called strong 
coupled regime, where crystallization or correlation can oc- 
cur. This regime as already been achieved for non-neutral 
purely ionic (or purely el ectronic) plasma in trapped system 
jPubin and 0'Neillll999() . 

The physics of ultracold plasma is very rich and 
only some aspects will be discussed here. Right af- 
ter the plasma creation, electrons which move faster 
than ions leave the sample within nanosecond time scale 



jTkachev and Yakovle'nkoll2001al) . This very complex stage 
has not yet been experimentally studied but is related to 
the Langmuir paradox which is merely the " violent relax- 
ation" process in astrophysics iChavani j|2002l) . This proba- 
bly leads to almost maxwellian distribution for electrons in a 
sub-nanoseconds time scale which is the in verse of the elec- 
tron Langmuir angular frequency ~ ^ 47r g|ne/ (eorrie). 
Similar coUisionless process occurs for ions, created in a spa- 
tially disordered state, but at hundreds of nanoseconds time 
scale giv en by the invers e of the Einstein an gu lar frequency 
loe ~ ^/A■Rqim/{eomi) l|Pohl et al.l l2004dO ISimien et all 
120041) . The laser ionization create a plasma from where elec- 
trons escape, creating a net positive space charge. When the 
potential depth becomes equal to the electron kinetic energy, 
occurring for an ion number equal to. 



N, = k2 
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1500 



(4) 



the space charge form a trapping potential for the new 
formed electrons. In fact the experimental results by 
iKilhan et a'i]|l999l taken for T;^ ranging from 4K to 800K, 
leads to a more complex result for the final p opulation dis- 
tribution. We found that lKillian et a'l]ll999l results can be 
well fitted (see figure |2J by the important formula : 



(-) 



(5) 



yielding to Ni ~ ~ 20000 < Ni which indicate that the 
plasma is usually quasineutral. We have obtain similar re- 
sults, but with a slightly bigger numerical factor in formula 
Equations Q s^nd can be written in a more conve- 
nient form: 



t:(k) = 



8.9 {N - Nef 
N 



(6) 



(T(^m) 

The electronic pressure finally lead s to expansion of 
the gaussian plasma. JKulin et alj l2000l have experimen- 
tally studied this expansion. They monitored evolution of 
the mean electron density fie using electron ejection in- 
duced by forced electron plasma oscillations (angular fre- 
quency tjpi = ^/ qifie / {some)) created by an external Radio - 
Frequency (RF) electric field iBergeson and SoenceillioO^) . 
The time dependence of the plasma size was then analyzed 
by Monte Carlo Method bv iRobicheaux and HansonI |2003 
and may be represented by a typical plasma expansion time 
tpE = ao/vo and an expansion given by: 



(J (t) ^ CTo + Vot 



2 
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ksT^ 
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Figure 1. Schematics of a typical experimental setup (typical 
size 10 cm) used to produced an ultra-cold plasma. The whole 
chamber is under vacuum only filled by the atomic vapor. The coil 
(radius of few centimeters), the trapping and repumping lasers are 
used to laser cooled the atoms. The pulsed dye laser excites the 
atoms into the ionization continuum or into Rydberg states. The 
electrons and (or) ions can be extracted by applying electric field 
on the grid surrounding the plasma. They are then collected by 
the charged particle detectors MicroChannel Plate (MCP). 




Ni/N* 



Figure 2. Experimental results bv lKillian et al Jl999l fitted with 
solid line by equation l5l . 



In fact after t ypically few microseconds an i onic den- 
sity spike appears (|R,obicheaux and Hansonll200,^ and for- 
mula ^ is no more valid but s eems to be restored for much 
long er time et al ]|20043). One major experimental re- 

sult l)Ktilin e^jT l200m concerns the fact that even for laser 
ionization at threshold, where a near zero electron velocity 
is naively expected, the velocity vo is no more given by for- 
mula (0 but is experimentally found to be always greater 
than 40 m/s and greater than yJUl jmi |Kulin ct al. 200^. 
One gets the fundamental result that the electron temper- 
ature Te is higher than 5K independently of the experi- 
mental initial parameters. The plasma is therefore always 
in a non coupled regime where < aws < ^d- This, can 
be due to instantaneous initial increase electron kinetic en - 
ergy while moving toward ions ( |Oizmii^Mid^]Neil||2002y) 
or du e to continuum lowering ( HahnI l2002l : iMazevet et alJ 
|2002^ which reflects the fact the energy of the isolated atoms 
is shifted by long-range Coulomb interaction with neighbors 



when embedded in a plasma. Another probable explanation 
results from huge three body recombinaison (TBR) rate. In 
such collision between two electrons and one ion, (Rydberg) 
atoms are formed in the plasma whe n one electron recom- 
bine with the ion. iKiUian et aT]l200ll have indeed observed 
in ultra-cold plasma this TBR. Rydberg atoms formation 
are the analog of binary system formation through three 
body encounters in globular cluster, this is a large energy 
source for t he free remaining electron whic h strongly heats 
the sample iRobicheaux and Hansoiill2002h . 

3.2 Three body recombination processes, binary 
system and Rydberg atoms 

The reverse process namely spontaneous evolution of an ul- 
tracold excited (Rydberg) gas to an ultra-cold plasma have 
been demonstrated by our group and by a group at Univer- 
sity o f Virginia simultaneously in year 2000 (iRobinson et alJ 
I2OOOI) . The Rydberg ionization process probably starts with 
blackbody photoionization or high energetic collisions with 
hot surrounding atoms. The so formed initial electrons 
leave, almost instantaneously, the interaction region. One 
attempted analogy with cluster might be dissociation of 
primary binary system due to external supernovae explo- 
sion but with much higher probability rate. A second phase 
occurs when the positive ion potential is deep enough to 
trap subsequent electrons, which then collide with Ryd- 
berg atoms cre ating more electrons in an avalanche ion- 
ization process iRobicheaux and Hansoi] l2003l : IPohl et al] 
I2OO3V This is perfectly identical to binary collision in stars 
clusters when collision with a third star leads to destruc- 
tion of the binary system. However, other relevant processes 
have been proposed whose effects need to be investigated , 
such as autoionization of Rydberg atom pairs jHahnll2000D 
similarly to destruction of binary s ystem through binary- 
binary collisions in globular clusters. IVanhaecke et al]l2004l 
have reported the total ionization of a cold Rydberg atomic 
sample embedded in an almost neutral ultracold plasma. 
This experiment is the analog of primordial binary system 
present initially inside a globular cluster. As in the astro- 
physical point of view this has strong influence on the ener- 
getic story of the sample. Rydberg atoms have finite lifetime 
(typically tens of microsecond) induced by photon sponta- 
neous emission or by blackbody photon absorption; these 
processes can be seen as the analog of stellar evolution. 
Some other works have shown that high-angular-momentum 
(circular electron orbit) Rydberg atoms, with order of mag- 
nitude longer lifetime (several millisecond), are created by 
eccentricities change th rough collision with electrons in an 
ultra-cold plasma (.Dutta et al.ll200ll : IWalz-Flannigan et alJ 
2004). This kind of process has been also mentioned in case 
of ZEKE (Zero Kinetic Energy) photoelectron spectroscopy 
(HQ ct al. 200 J). The high-angular-momentum Rydberg 
atoms are the analog of circular bin aries in cluster envi- 
ronme nts which is still an open st udy (iGiersz and SpurzemI 
20041) . IWalz-Flannigan et al.lf2004l also note that there is an 
unexpected large population of deeply bound (hard) Ry- 
dberg atoms, this is probably due to Penning ionization 
where two binaries collide leading to the disruption of one 
of them and increasing of the binding energy of the second. 
This Penning ionization is the exact analog of the disrup- 
tive collision in binary-binary scattering in the cluster case 
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jLvman SpitzeillQSTtlHut et alJl993) . Indeed, the coUisions 
during binary-binary and binary-single interac tions still 
need a lot of investigations ijFregeau et alJl2004h . Another 
interesting study is relative to the binding energy distribu- 
tion. Starting from a pure plasma the binding energy distri- 
bution, resulting o f the TBR process, is modified during the 
plasma expansion iRobicheaux and Hansonl2002f) . The cen- 
ter of mass Rydberg velocity is the one of the ion at the time 
the recombination occurs and should therefore be linked 
with the Rydberg binding energy : the slower atoms tend 
to ha ve larger binding energies iRobicheaux and Hansoiil 
120031). This is also linked to the so-c alled excitation freez- 
ing jTkachev and YakovlenkdEoOllJ) and needs to be ex- 
perimentally studied. There is however di fferences between 
Rydb erg atoms and binary stellar systems llHut and Bahcahl 
Il983r (see also the nice revi ew on binaries in globular stars 
clusters bv iHut et al.l[l992l) . Indeed, Rydberg are dipolar 
but neutral atoms. Therefore the Rydberg binary system in- 
teract with neighbors only through charge-dipole or dipole- 
dipole interaction scaling respectively as and a^g but 
in the cluster case binary system interact through Coulomb 
interaction scaling a^g. Furthermore, the dissipation of 
the relative kinetic energy of two strongly interacting stars 
through tidal effect which lead to the formation of tightly 
bound binaries (so called tidal binaries) has no simple ana- 
log in the plasma case. The goal of the article is to give an 
overview of the analogy between both systems, but it is be- 
yond its scope to detail all the process involving binaries. 
We will mainly focus on results which are of interest for the 
free electrons distribution. 

When a Rydberg gas is formed, there is competition be- 
tween the deexcitation rate and the excitation rate for Ry- 
dberg atoms. In fact, very highly excited Rydberg atoms, 
for which binding energy are less than the single electron 
plasma kinetic energy, are analog to soft binary stars and, 
as in the cluster case, encounters lead to disruption of 
the system by gradually increasing the binding energy un- 
til it becomes positives. This is a similar law than Reg- 
gie's law (we will use this terminology even if this law 
was known in atomic physics much earlier than in clus- 
ter physics ) : hard bina ries get harder and soft binaries 
get softer iHeggid Il975^ . Indeed, the probability of exci- 
tation for a Rydberg atom equals the probability of de- 
excitation when its binding energy is given by 3.82fc_BTe 
(Mansbach and Keck 1969; Robicheau x and Hanso n 2 0d3) . 
A similar formula but with a bottelneck point nearer t o 
fcsTe h as been given b y |Tk^hgv_aml_Vakovler^^ (|2001b|), 
see also IStevefeh et al.l ^1975^ or Tvri^n^m^^m^etj lll98(J) 
for other formulas. This is exactly the same behavior than 
the one noticed in stars cluste r physics. This pro c ess have 
been studied experimentally iKillian et alJ I2OOII: iLi et all 
l2004l:IWalz-Flannigan et al.ll2004) . 

As previously mentioned, in both systems, the binaries 
formation comes from thr ee body recombin aison (TBR) . For 
clust er, the rate Ftbr is iBinnev and Tremaineiil987.) (8-7) 
and llLvman SpitzeI^ll987^ f6~3r) 



where = ksTe/me is the one dimensional ve- 

locity dispersion. And for electrons plus ions we 



have iMansbach and Kecklll969l : iTkachev and Yakovlenkol 
l200ld) 



1 / drie 
rie \ dt 



~ —Ftbr 
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In both formula we have omitted numerical factor of the 
order of unity. The analogy is obvious but the T^^^^ be- 
havior leads to a main difference between the cluster case 
and the ultra-cold plasma case. For an ultracold plasma 
the TBR rate is huge but it is almost negligible for clus- 
ters. We mention that the analogy could be pursuit further 
on. Indee d, formula given by (Sigu rd sson and Phinnc y 199 J 
(see also iHut and BahcalJ lll983t ) : TCvman Spitzeif 11987)) 
which concern the energy change between sta rs is similar 
to the rate function jMansbach and Kecklll969h (in-12) be- 
tween atomic energy states £i — kBT^ei to energy state 
£f — ksTeef given by: 



k{€i,e;) = 

k{ei,€f) = 



fco(-e,)'"^(-.,) = 
fco(-eO-''^^e-(^'- 



(9) 
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where fco = 11 y ^.^^ j ( In" ) ^'^'^ energy is 

taken as the ionization threshold. 

In the cluster case, the energy heating rate due to bi- 
naries formation is 



£ 



100rTBR?n(T„ 



the factor hundred (which is a very approximate one) comes 
from the fact that every binary liberates in the cluster 
an energy of few hundreds times ksT through encoun- 
ters with other stars before being ejected ou t of the clus- 
ter by a very energetic reaction jCohn et al.lll989ll . This 
non-negligible fiux of escaping binaries is a we ll known ef- 
fect in cluster case iMevlan and Heggid Il99'^ . But, flux 
of escaping binaries is probably negligible in the plasma 
case, because the binding energy of the binary increases 
mainly due to radiative lifetime and not much due to col- 
lisions. The bottleneck occurs when its energy becomes 
lower than e* » fca x 500 K x (-i%l 



-2/9 



,1/9 



J ~ II iviiv/ V 10^ cm / 

iTkachev and Yakovlenkdl2001br) . Therefore in the plasma 
case, we do not have such simple results for £ because the 
radiative lifetime combined with t he expansion play a com- 
plex role, but we could estimate iTkachev and Yakovlenkol 
[2001b) (12): 



£ ^ 5.4 



109cm-3 3/is 



-2/9 



TBRfcsTe 



(10) 



The analogy holds because the process evolve as £ oc 
rxBRfesT in both cases. 

Finally, it might be useful, ofr further studied, to 
note that if the time dependence Rydberg bi n ding en - 
ergy distribution is complex (see IPohl et all (|2004bll : 
Robicheaux and Hanson ( 20o3)), we found that the results 
can be well fitted by a binding energy ^Ryd > distribution 
of the Rydberg states proportional to 



-fRyd/C^flTRyd) f '^Ryd 

fcsTRyd 



(11) 



with a and TRyd are time dependent parameters and a is 
always smaller than its equilibrium value of 5/2. 
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4 FOKKER-PLANCK EQUATION 
4.1 Symmetry consideration 

iRosenbluth et al]ll957l already indicated that the Fokker- 
Planck equation is formally identical for electrons in plasma 
or for stars in stellar systems. We thus expect similar be- 
havior for our system or for an isolated cluster. Velocity 
anisotropy is generated by t wo-body relaxation in the outer 
part of a globular system jLvman Spitzeil Il987^ . We ex- 
pect similar results here but, to simplify the discussion, we 
will use isotropic symmetry. Using Lynden-Bell's improved 
strong Jean's theorem for a spheri cally symmetric plasma 
system jBinnev and Tremainel[^987^ or assum ing ergodicity 
in ou r ionic non harmonic trapping potential iSurkov et alJ 
Il996l) . we could assume further on that the electronic phase 
density function f{r, v, t) in the ultra-cold plasma depends 
only on energy £ = qE4>+-^meVe where cf> = (f>e+4'i is the sum 
of the electronic and ionic potential. To be able to compare 
with stellar dynamics it is more convenient to use energy per 
mass notation E — S/nie — $ + ^v^. The potential energy 
(per mass units) is $ = "3>i -I- $e = ge(0i + 4'e)/'me- Using 
the Poisson's equation one gets: 

Or \ dr J eo 
N,{r)-N,{r 



A* 
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dr 



dr ) 



A-KG'(pe+Pi) 



(electrons or ions). Tf — 47rG^mjlnA for clusters becomes 

2 2 

Ff — 47rG''^m?lnA — . 'I' ., InA = F in the plasma case 

which, in this case, is independent of the field particle mass. 

If we use the so called thermal bath (i.e. Maxwellian) 
approximation for the field distribution function //, one 
gets: 



dt 



E 



dv 



G{ 



) f 



ruf 



+ 



(vj) V df 
3 dv 



(16) 



Where G{x) = Erf(a:) 



-(s) tends to 1 when x tends 



> 



pi 
(12) 



to infinity and is proportional to x for small x values. Us- 
ing relations Ti < Te (or {v}) < (wf)) and < rie 
we found that in equation (I16II the coUisional ionic terms 
are negligible, in the field particles coUisional part, com- 
pared to the electronical ones, with Fokker-Planck of pure 
electrons. CoUisional ions effect can be summarized mainly 
as added to ne and a {vf) terms added to = vi 
in the electronical part. The arguments hold also for non 
maxwellian distributions and the key results is that we can 
always safety neglect ions in the coUisional part. Finally 
the (L andau)-Poisson-Fo kkcr-Pla nck equation for electrons 
is then jPelcroix and Bers 1994: Mitchner and Krugel^l992^ 



9/ 
dt 



dl 
df 



df 



dl 

dv 



47rr d 

Il2 Qy 



V^ 



f f'v"dv'+ 



$(r) = $oo - 47rG' 



(r )r 







(r )r dr 



where pe = men, 
Ps + Pi , and Pi — 
Pi = miUi 



is the mass density for electrons, p° = 
-meUi which is not the ionic mass density 
The artificial "total" mass is M*°* = (Me + 



2i^M0 = m^Ne - N,) = + M, < 0. iVe(r) = M./m, 
(respectively Ni{r)) is the number of electrons (respectively 
ions) inside a sphere of r radius: 



vdf ( r w% , 



f'vv'dv' 



(17) 



the prime indicates that the quantity depends on the en- 
counter field velocity v' , not on v. This equation contains 
only electrons and is thus completely identical (using equa- 
tion Q) to the Fokker-Planck equation for a single mass 
stars system. 




N,{r) -- 

Erf is the error function and 
Erf I 



<E>i(r) = Gm^Ni 



(13) 



(14) 



V2rT 



The magnetic field (see coils in figure is usually 
turned off during an experiment or its effect is small and 
can be neglected in first approximation. The Fokker-Planck 
equation for the electron space phase density distribution 
f ~ fe, which can be seen as a series in In A and is then 
no more valid in the strongly coupled case (A < 1), is 
jPelcroix and Bers 1994; Mitchner and Kruge3 ri992ll 



dt ^-^ v^ dv 

f 

vdi 

3 dv 

with 



ffA-KVjdvf -\- 



(15) 



v'f ff^nVfdvf + V 



-^ffi-Kvfdvf 



^^21^^ di_d^ dl 

dt dt df df dv 
where the subscript / indicates the type of field particles 



4.2 Virial theorem 



There is fundamental uncertainties in the value of In A. For- 
mula InA = ln{8nq~'-'{£okBTef^^{n'i] 
for homogeneous systems wher e ne(r) 
tem an usual approximation is llLvman SpitzeJl98i1) (2-14) 
A ~ 2 ^'^^'^ where rn the radius containing half mass of the 



^^^) is only correct 
n^. For a s tar sys- 



system and f_L = . Similarly for plasma we define a 
(global) temperature Te jBinnev and TremainellT987ll by: 



2Ke = 3kBTeNe 



d'^fd''v — meiff{f, V, t) 



where Ke is the total electron kinetic energy. By anal- 
ogy, in the plasma case, we can write A = 2^ with 
= —G'ml/kBTe. The equations are then perfectly identi- 
cal for both systems except the presence of the Debye screen- 
ing distance in the Coulomb logarithm for plasma but the 
radius limit distance for cluster. One usual average value, 
used in single mass globular clusters, consists to evaluate 
using the potential energy 



tot J- 

W = — 



47rrV*°*(r)$(r)dr « -0.4-:^ (18) 

rh 



where the constant 0.4 is a reasonable approximation for 
most systems (0.44 for gaussian density distribution where 
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Th ~ 1.54cr) iLvman Spitzejll987^ . The static Virial theo- 
rem 2K = — leads to A « QAN where is the number 
of stars in the cluster. In the plasma case the Virial theorem 
iBinnev and Tremain e 1987) is: 



2K, 



pe{r)r 



dr 



■Anr dr 



We: 



(19) 



Where Wei = /q°° Pe(f )f ■^|r47rr^cir can be seen as an ex- 
ternal "potential" energy (due to the ions) and Wee ~ 
I J pe$e. A naive extrapolation, based on quasineutrality 
approximation (n^ — Ui constant), of the cluster case would 
leads, in the plasma case to: 



2Ke = 3kBTeNe 
Te{K) 



0.4 



(20) 
(21) 



which have to be compared to formulas 1181 . (|SJ. Finally, 
using this naive evaluation of the Virial theorem we could 
extend for the plasma the formula A ~ 0.4A'^ valid for clus- 
ters in: 



A = o.4(Af, -iV,)^ 



(22) 



Q.1{N, - Ne) 



Ni-Ne 



Ne 



where we have used formula 12111 and the gaussian approxi- 
mation Xd ~ \/ {2-Ka'^)'-^/'^eokBTe/{qlNe) to derive the final 
formula. We will see later that electrons are not in a gaus- 
sian distribution and that these naive formula have to be 
corrected. Nevertheless, they indicate three very important 
results. Firstly the ions minus electrons number is directly 
related to the temperature. Secondly, the final temperature 
Te is on the same order to the one expected from the laser 
wavelength T2 . Thirdly, the coupled regime (A < 1) can be 
reached only for almost pure neutral plasma with numerous 
electrons. 



5 FLUID AND GAS MODELS 
5.1 Gaseous equations 

One of the easiest way toward an approximate solution of 
the electron Fokker-Planck equation 11711 is to use the veloc- 
ity moment equation. Here again, to avoid confusion and to 
be able to use equations derived for stars dynamics we define 
the mass density distribution fm = nief ■ It is simple to use 
dimensionless equations where r = vor* ,t = tot* , p{r,t) = 
pop*{r*,t*),Mr{r,t) = MoAf*(r*,t*), . . . or the selfsim- 
ilar form where r — re(t)rt, p(r,t) = pc{t)p,{r,), . . . 
iLouis and Spurzemll99ih . We then define the mass density 

p = Pe = rUeUe = (0) = / (fvfm{v) = °° A'Kv'^ fm{v)dv 

related to Mr = meNe{r) the mass contained in a sphere 
of radius r. We also define the velocity of mass transport 
u hy pu = (1) = J Vrfm{v)d^v where Vr is the radial 
velocity, the kinetic energy density ks by 2kE = (2) = 
°° A'jTV^v'^ fm{v)dv. This also define the pressure p by 
3p -|- pu^ = 2kE which is linked to the temperature Te, the 
one dimensional velocity dispersion a„ (isothermal sound 
speed) through p — UeksTe = pa^. 



The moments equations are the following: 
• The mass integration 
dMr 



dr 



= 4-Kr p 



which can be written as 

91nAf* .2 , 

ir-\ — T = ^ P 
omr* 



for the choice Mo = A-Kr^po 

• The continuity equation (with mq = fo/to) 

dp , , dp 1 dr^pu 

= 7^+div(pn) = -^ + 



dt 



dt 



dr 



(23) 



(24) 



(25) 



dp* u* dlnp* du* /r* „u* 

= 7577 + — TTT-^ + TTT^+S- 
OT r* omr* dmr* r* 

• The hydrodynamical (Euler) equation 

Du 1 dp 9$ 
Vt p dr dr 

where = ^ + is the co-moving (Lagrangian or con- 
vective) derivative following the mass evolution ^^^f- ~ 0. 
The third momentum equation is : 

• The kinetic energy transport equation 



(26) 



dks , 1 , 9$ „ 

— -f-div{3) + p^.— =0 



(27) 



The fact all these equations do not contain right hand 
side collisional terms is a verification that the Fokker-Planck 
equation conserve the mass and the energy. Therefore there 
is no difference between the use of the collisionless Boltz- 
mann (also called Vlasov) equation and the use of the colli- 
sional Fokker-Planck equation at this stage. 



5.2 Self similar collisionless quasineutral plasma 
evolution 

The main goal of this paper is electrons evolution, but simi- 
lar gaseous equations holds for the ions in the plasma. In the 
ion case the pressure is negligible {Ti ~ 0) and the system 
of gaseous equation is closed with 



Dui ^ 1 dpi ^ Qi d(j> ^ Dui ^ Qi me d$ 
Dt miUi dr rui dr Dt Qe nii dr 



The quasineutrality equations: qiUi + qeUe ~ and Qi-^ + 
qe — can be used, with the Euler equations for ions and 
electrons, to lead to a useful relation between the potential 
and the pressure : 



dr 



dr me dr 
qln^ _|_ "^f"' 



(29) 



Assuming Maxwell-Boltzmann electron distribution we 
found 

_^(,,t)^,M£W 
dr ' mia{t) 

The dimensionless equations are ver y similar to the self 
similar ones used in clusters physics l)Louis and SpurzemI 
1991). The selfsimilar ionic evolution is easy to determined 
and verify: iRobicheaux and Hansonll2002L^2003^ : 
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Ui{r,t) ^ tr- 



2 2 



1 + 



(30) 



tion (Tfll . iDorozhkina and SemenovlllQOa have demonstrated 
that this result is in fact more general. Using Vlasov equa- 
tions for /e and /i, assuming a quadratic form for the po- 
tential, a self similar homological evolution, quasineutrality 
and equation l|29|l they also found equation l|8()^ . In fact 
this result has been recently improved and the most gen- 
eral form for a self similar quasineutral coUisionless plasma 
expansion into vacuum, witho ut any other assumption is: 
iKovalev and Bvchenko'vll2003h : 



+ 



is 



(31) 




r/o(t=0) 



Figure 3. Ion density evolution based on equation I32i for u 
250Atm and Ne = 20000. 



Using n = ^ we recover equations l|3()|l and Q which are 
then shown to be the most general self similar solutions for 
a quasineutral plasma. 



5.3 Temperature evolution 

One important question concerns the temperature evolution 
in this self similar solution. The energy conservation comes 
from the time derivative of the square of equation I12II (Am- 
pere's law): 

__._(p,,,+p,„0 
added to equation I27I I and spatially integrated: 

-ksTeiO) = -kBT,(t) + -m, J u'^,(r,t)n,{r,t)4nr'^dr 



where we did not take into account the negligible ionic ki- 
netic energy Tj <C neither the ionic correlation energy 
JPohl et al.ll20 04df^ or the heating due to the three body re- 
combinaison (see Robichcaux and Hanson (2003) and equa- 
tion lion '). We see that during the plasma expansion an elec- 
tronic adiabatic cooling occurs and the ionic distribution 
function becomes a Schwarzchild-Boltzmann one given by 
equation 13111 . 



5.4 Analytic ion spike evolution 

The coUisionless theory, we have just developed, was based 
on several assumptions and has the usual drawback of all 
self similar solutions predicting unphysical results as, for 
instance, a velocity increasing without limit for r going to 
infinity. Furthermore the theory is only valid before the ionic 
spike appears. 

We know from Monte Carlo ijRobicheaux and HansonI 
l2003 h and molecular llPohl et alJ2004b^ simulations that the 
quasineutrality is violated and that the ionic front shows 
a density spike after few m icroseconds. Some theoretical 
prediction (|Pohl et al.ll2004'b|) predicts that the ionic spike, 
where the quasineutrality is violated, disappear. This result 
has to be related with the work on non-ultracold plasma 



where it can be shown that the ion front m oved at velocity 
varying loga rithmically in time llMorall2003D . 

FoUowing iKaolan et ai]l2003l. who studied a non neutral 
plasma, we will describe the r-radius (with Ni{r) ions) ion 
shell evolution. We assume that along a r-shell trajectory 
r{ri,t) starting at the initial point r^, the total number of 
ions remains unchanged: Ni{r{ri, t)) — Ni{ri). The ion den- 

sity is then given by (r, t ) = ni{ri,t — 0) . The condi- 
tion that no particle trajectories cross each other is violated 
when the shock formation appears creating an infinite lo- 
cal ionic density. A natural expansion of the quasineutrality 
is iVe(r) « Ni{r), we will therefore assume Ne{r{ri,t)) = 
Ne{ri). The ion Newton equation is: 



47r£o 



integrating ^ times this equation leads to the conservation 
of energy equation and to the analytical implicit solution: 

t 

^3/2 y "47reo 



(32) 



A typical result with gaussian a ~ 250/im approximation 
for rii and Ue and with Ni — Ne = 20000 is given in figure |3 
The typical "Coulomb Explosion" time scale is then 

~ 3.5 fis where — n{ri — 0,t = 0). 




tcE 



g'i 



At early time t <^ tcE the evolution equation 1132^ leads to 
an evolution given by 



- = i + r^-A/- 



4 V 2 



where we have used equation With equation and 
with the gaussian approximation 113|l this might be written 
(for r — a): 

0.8(f)"%'| ,33, 



-«1 

O"0 



which is, with (^)°^ on the order of 10, is very simi- 
lar to equation (|7| for early times. Even if the hypothesis 
Ne{r{ri,t)) — Ne{ri) is probably worse than the one for 
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ions. The final results should not be too badly affected by 
this assumption. Indeed, the only requirement for equation 
is the time independence of Ni{a{t)) — Ne{o'{t)) which 
is probably a better assumption. 



5.5 Closure relations for the electrons 

Let us go back to the electron evolution to find the time 
dependent evolution equations for the electron distribu- 
tion. According to the ambipolar diffusion and due to the 
quasineutral plasma behavior we have (especially in the core 
region) u ^ m so (in absolute value) 

Du ^ Dui 
In ~ ~Di" 



defined for local homogeneous system. Using dimensionless 
values we can see that to — teo = 9a„o/(16v'7rG'^mep° In A) 



and I, 







1 2 



3o-„ 



3^2 



This leads to formula: 



Ac oc — = 



--^3(-G')me— InA 



(38) 



me 9<1> 1 dp 
rrii dr p dr 



where the dimensionless proportionalit y factor shou ld be on 
the order unity (nearly 0.4 = g^C" in lLouij il99Ct) '). 

The phenomenological heat fiux formula (13 711 is based 
on the assumption I > a which is not always verified, espe- 
cially for small Coulomb logarithm value. Furthermore the 
fact we use only local variable for the thermal conductivity 
equation (with vq = ro/to): 



Therefore we need to go at least one step further to close 
the moment equations for the electrons. For our isotropic 
quasistationary evolution 3p = (2) and we define (see 
iLouis and SpurzemI llGQlh ) the velocity energy transport 
V by 5pv = (3) = J Vrv'^ fm{v)d^v which is linked to the 
heat flux (different from the energy flux across a constant 
radius sphere) F — 5p{v — u)/3 = L/(47rr^) where L is 
the luminosity. The hydrodynamical equation becomes with 
M* = m^Ni/Mo: 



V u 

r* r* 5A 



9 p* 



^i/2p*i/2 a(inp* -inp-) 







^1 * " (39) 

migh t not be the best solution (see iLouis and SpurzemI 
119911) 1. It is also po ssible to use t he hig her momentu m equa- 
tions, calculated bv lLarsonll970l and bv lLouiJl99'ol namely: 

dr ^ dr 

where 15k = (4) = J(^°° 'iivv'^ fm{v)dv and 



(40) 



dlnr* p* r* 



(34) 



- + -div(.(2«-u)) + -p- 



5U 



(n-p'/p) (41) 



where A — Ne/{Ni — Ne) and we have assumed po — 
{—G'){Ni — Ne)mePo/ro ~ poOvo^ ■ A = —1 is the case with- 
out any ions and, in this case, equations are identical to the 
cluster case. 

The energy transport equation becomes 

dp 5 , , 2 9$ 

91np* 5u* 91np* 5 dv* /r* v* 2 u* p* M* 

= + ^, + TT^TT-^ h5 -A- ^ 



dt* 3r* 91nr* 3 91nr* ' " r* 3 p* 
which can be written in a first thermodynamical law form 



We have shown that the equations are almost iden- 
tical in cluster and plasma systems with an extra A fac- 
tor in the plasma case. We shall not fully solve the equa- 
tion here but the classical iHenvev et alj il964l) -Newton- 
Raphson implicit difference relaxation method is probably 
well adapted. In fact we shall not resolve here the gaseous 
(^ggl^ation because the original Fokker-Planck equation should 
give better results. Let us note t hat using A bel's transform 
' Wp(<i>(r)) = 2^5/^ f^{E)yjE -^{r)dE it is possible 
to recover the phase-space density distribution function: 



dL , 2 \ D SkeT, , D 1 



(36) 



ME) = 



d 



VSn^ dE 



dp/d<t 
V<S>~ E 



The most delicate assumption is the form of the clo- 
sure relation. A simple one is the thermal conductiv- 
ity assumption given by jLvnden-Bell and EggletonI [T980t 
iLouis and Spurzemlll99lll : 

• The thermal conductivity closure relation 

F = ^K^^-A^^ (37) 
dr dr 

where K = kBAc/rUe is the thermal conductivity coefficient. 
The theory of heat fiux in gases indicates that Ac should be 
on the order of pl^ /t^oi- Where I is the mean free path and 
tco\ is the time between collisions better to be taken as the re- 
laxation time te (with velocity a-u ) see jLvman SpitzeJl987l: 
|L^id«i^Bel^n^ EggletonI ITgSoll for the cluster case and 
jMaiumdar et alJll973l'r for the plasma case. When A is not 
too small, the mean free path I ~ Stefi, (typically 2 mm) is 
comparable or larger than the size of the sample a. This indi- 
cates that the electrons make severals orbits before colliding 
with neighbors. Hence, it is better to use for I a typical ra- 
dial distance Ir between encounters given for in stance when 
the mean velocity (t„ is the circular velocity Ir 47r(— G')p/3 



This is the well known Eddington formula (1916) 
iBinnev and Tremainelfl987^ . 

The simplest solution we might think is a self similar 
one for electrons as well as for ions which in fact lead to a 
stationary solution. Gravitating systems never form static 
homogeneous equilibrium. However, plasmas contain both 
positive and negative charges, so they can form static equi- 
librium. The time scale for ion expansion is order of mag- 
nitude the electrons relaxation time toward equilibrium. So 
electron will reach equilibrium on a nanosecond time scale 
where the ionic motion could be considered as frozen. We 
might then first look for static equilibrium distributions. Un- 
fortunately equations 1361 leads to L constant and equations 
1371 1 and 1381 to L{r = 0) = and (t„ constant. A stationary 
solution has no heat fiux and is isothermal with a density 
distribution given by: 



d_ 

dr 



pi) 



This leads to a density distribution with a r~ power law at 
large distance meaning an infinite mass. In fact such distri- 
bution present an attractive potential for electrons up to a 
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point rt where the ion and electron density becomes equals 
and electrons evaporates leading to a repulsive potential for 
r > rt. We a lready know from Monte Carlo simulations 
iRobicheaux and Hanson 2003) that this isothermal solu- 
tion is correct for the core system but the finite mass "halo" 
can only tend to be isothermal. 

Before going back to the Fokker-Planck equation, we 
can mention that it is possible to take into account the three- 
body collisions through the heating term (IIUI in the kinetic 
energy transport equation. We should also use equations @ 
and IIH to give a full picture of free and bound electron 
distribution. 



6 ORBIT AVERAGE OF THE 

FOKKER-PLANCK EQUATION 

6.1 Truncated assumption 

If the system is isolated the maximum energy above which 
one the particles are extracted from the sample is Et — 
$(r- = oo). But, the system can be truncated at a given 
(apocenter) distance rt and the truncated energy is then 
i?t = $(r = rj). In a isotropic Fokker-Planck equation one 
has no other choice but to use the energy as a criterion for 
escape. For cluster physics the truncation in energy is a good 
approximation even if it is known not to be a sufBcient cond i- 
tion to describe accurately the physics jKim and OhllTooi^ . 
For cluster the differential forces, produced by the galaxy, 
result in two saddle points (distance from the cluster center 
which is then proportional to the cube of the mass) through 
which stars can pass over into the galactical field. For plasma 
physics under magnetic field a velocity truncation is proba- 
bly another good choice due to the velocity dependent force. 
But, the experiment can use an homogeneous electric field 
(see figure Fcy along the y axis. In this case the trunca- 
tion arise only from a single saddle point at yo = rt abscissa 
along the y axis where 

F = -^(yo) = --^(yo) (42) 

Therefore the truncation is not at all spherically symmet- 
ric and the hypothesis is probably worse for plasma than 
for tidally limited cluster. Nevertheless, under ergodic as- 
sumption, the hypothesis of truncation in energy might be 
a good one and will be used hereafter. There is always a 
stray electric or magnetic field in the experiment. The sys- 
tem is therefore never perfectly isolated and cannot extend 
to infinity without touching some electric grids (see figure 
Ql. Furthermore, an external electric field F is sometimes 
constantly applied (as done bv lKulin et aljIioOOt) to be able 
to accelerate electrons (or ions) toward a detector and to be 
able to monitor the behavior of the sample. In such experi- 
mental conditions equations 1421 . I12II and 11311 imply 



F 



g"- 




G'm, 



' t 47r£o ' t 

This experimental electric field can be experimentally 
used to extract electron trapped in an ionic cloud in order 



(43) 



to simulate the tidal escape from globular cluster under ex- 
ternal galaxy gravitational attraction. Typical value, before 
the expansion, are rt ~ 1.5 mm« 6(j for _F = 1 V/m. We can 
see here that rt is proportional to the square of the mass and 
not to the cube of the mass a s in the star cluster c ase under 
tidal effect of the host galaxy iLvman SpitzeIll987^ . Assum- 
ing the system reacts on the external field we will neglect its 
role inside the system simply summarize it by the Et value. 
This assumption is better for a small external field. In our 
experiment the magnetic field gradient [B' = 1.5mT/cm) 
is not turned off and it ejects electrons for r > rt = 12a 
(value estimated for Te ~ 100 K). This value is similar (see 
equation 1431 to the radius found for electric field on the 
order of _F = 2mV/cm which is indeed a typical value for 
uncontrolled straight electric field. As already discussed, the 
mean free path is usually large compared to the sample size, 
or similarly the orbital timescale at rt is negligible com- 
pared with the relaxation time. Thus, we will assume that 
the electrons with energy E > Et would be lost almost in- 
stantaneously. We k now this criterion is strongly violated in 
the star cluster case jTakahashi and Zwartlll99^ but this is 
nevertheless a reasonable choice, yielding to f{E) = for 
E > Et and by continuity 



f{Et) = 



(44) 



The zero probability presence at the border is oppo- 
site to the conclusion obtain in ultra-cold atomic system 
interacting through van der Walls atomic interaction po- 
tential. Here the Coulomb (or Newton) interaction acts at 
very long distance and encounters are "gentle" leading to 
a small continuous variation in energy during collisions. On 
the contrary the van der Walls interaction is at so short 
range that only the strong collisions with a single encoun- 
ters are involved. Consequently, in the route toward quan- 
tum degenerate dilute ultra-cold atomic Bose Einstein Con- 
densates the evaporation (it is more an ejection process but 
this "atomic" terminology is well known) process leads to 
truncated maxwellian distribution with non zero probabil- 
ity presence at the truncated border. Despites the differences 
similar strategies are used in stars system and in atomic sys- 
tem : superstar strategies in cluster code and macro-atoms 
strategies used in evaporative cooling models of ultra-cold 
neutral atomic cloud in route to Bose-E i nstein condensa tion 
jBerg-Sorensenlll997l: iPinke et aljll998l: iTol et aljliooil . 

6.2 Orbit averaging 



Using V = -y/ 2{E — "!>) we define the volume r of phase space 
with energy less than E and the isotropic average q of the 
radial action by: 

q{E,t) = i /" ' \2(E-<f{r,t))'/^r^dr = ^r(£,t)(45) 

For our quasi gaussian system we found, from equation 1141 
that a good approximation (within 10% accuracy) for q is: 

\3 



q{E) oc {Et - E){E - EoY 



(46) 



This is not surprising because we are between a square po- 
tential, where $ = and q oc E^^^, and a quadratic poten- 
tial, where $ = $(0) + iojV^ and r = ^{E - Eo)^. We 

dq 



will also sometimes use 



/; ^"VV2(i5-$(r,t))dr 
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which is the isotropic average of the period and is 1/{16tv^) 
times the classical microcanonical phase space density func- 
tion. Indeed, on the contrary to the ultra-cold atomic sys- 
tem case, the plasma is far from quantum degeneracy - the 
Fermi temperature is on the order of tens of microKelvin 
which is orders of magnitude below the plasma temperature 
- so we can use the classical density phase distribution func- 
tion which is the same than the quantum one but multiplied 
by the fundamental phase volume {h is the Planck's con- 
stant). 

Multiplying the Fokker-Planck equation 1171 by the 
delta function S{E — + $)) and integrating over the 

phas e space leads to the orbital average Fokker-Planck equa- 
tion jHeno Jl96 J: iLvman Spitzeilll98 11) 



dE dt dt dE dE 



\dt)r 



21 

dE 



f'q'dE' + q / f'dE' 

En J E 



(47) 



(48) 



where Eo = $(0, t ). Equation II48[I is merely equ ation H47|l 
divided by dr/dE ilnaeaki and Lvnden-Belil99Cl) . 



[). The flux 

through phase space 11 can be written in another nice form: 

nEt 

n = 47rr / 

J En 



"//'Min(.,.') ('J^^'J^)dE' 
^ ' ' \ dE dE' ) 



which deflne the generalized concept of temperature T^{E) 
by 

1 dlnf _ 1 

~me dE ^ kBTp{E) 

The formulas can be checked with the use of a square poten- 
tial ($ = for r < rt). We have, in this case as in free space, 
q oc E^^'^, E = and equation H47|l recover exactly the 

Fokker-Planck equation H17^ . 

Equation 14711 is called the orbit average Fokker-Planck 
equation and holds when the relaxation time is longer than 
the crossing time or the orbital time. As previously men- 
tioned we can experimentally violate this condition (for very 
cold electrons or very large and dense sample) but usually 
this assumption holds. Thus the distribution function will 
evolve slowly compared to the orbital period of each elec- 
tron. In this slowly time variable potential (orbit-averaged 
assumption) 



dE/dt : 



drrv{d^/dt)/ 



dq/dt 
dq/dE 



hence 



and 



adia- 



dt ' 

^ = and therefo re f(E,t) = 

iBinnev and TremaindligsV : Inagaki and Lvnden-Belil 



is the mean value of 
batic invariant: 4t 



Methods of calculation, developed for stars clusters, 
can then be very easily adapted. For instance the Chang 
and Cooper (1970) finite-differencing scheme used in Cohn's 
method (Cohn 197!^ is perfectly adapted. Indeed, the first 
step of the method is to advance / in time $ being held 
fixed. The Fokker-Planck equations are identical for both 
systems so no change has to be done in this code step going 
from stars system to ultra-cold plasma system. The second. 



easier step, is to advance $ by solving the Poisson's equa- 
tion with / being fixed as a function of the phase volume 
adiabatic invariant r. Other two step methods can also be 
used, as the one developed by Jakahashi 1993 from varia- 
tional principle. We hope that the simplicity of this project 
would lead to rapid use of Fokker-Planck code in ultracold 
plasmas physics. 

6.3 Evaporation and ejection 

Electrons can escape from the plasma by ejection (through 
strong collisions with a single encounters) or be evap- 
oration (through series of weaker distant encounters) 
iBinnev and Tremainelfl987l) . 

The eva poration rate has a long history in stars dy- 
namics (e.g. llJohnstonelll993Dl Let us first formulate, for 
our plasma system, the THfaonl (jl960) paradox known in 
clusters physics: in the case of isolated system (rt = oo) 
an electron with energy E ~ Et spend most of its time far 
from the center and suffer very few encounters. Its change 
energy rate goes to zero when E approach Et so the electron 
never escape in this diffusion type (Fokker-Planck) picture. 
Thus, the orbit average Fokker-Planck equation, with the 
same relaxation time f or all electron s, do not hold for outer- 
most "halo" electrons iSpitzer and Shapiro 1972). The fact 
that halo electrons play a role has been in fact experimen- 
tally sk etched using R adio Frequency heating of electrons. 
Indeed. llTet alj|2004l used RF electric field to "shake" the 
trapped electrons and speed up the ionization of Rydberg 
atoms inside the plasma. They have noticed that the pro- 
cess is less efficient for large RF amplitude due to the fact 
the forced RF electron oscillations drive electrons outside 
the cloud from the center reducing the amount of time the 
electron could collide with Rydberg atoms. This experiment 
indicates that when electrons are in the halo they do not 
collide that much. 

From Henon's paradox we can say that the electron can 
escape from an isolated system only from ejection due to 
strong single close encounters. The rate has been calculated 
by Hcrion 1960 in the cluster case. By analogy using equation 
Q we could then assume: 



dN^ 2, , 2 2n2 

— = -(16Gm..r) 



r drx 



{Et - E)^ J\ )J\ I 



(49) 



where the integration hold for E' ^ Eo, Et E Eo and 
E + E' ^ $ + Et- Performing the integration over E' leads 
to the loss rate of electrons with energy E. 

The case of a truncated system {rt < oo) is different 
and the evaporation process is a diffusion in velocity space 
process well described by the Fokker-Planck equation. The 
number N{E) of electrons with energy less than E is 

N{E,t)= f f{E',t)^{E',t)dE' (50) 



time derivative straightforwardly leads to 



When E - Et {or oo),N - and the first (atrd third) 
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"spilling" terms are null (see equation (1441 ) . The second one 
can easily be calculated using the general form 

fiE,t) = f(E,t)Y{Et{t)-E) 

where Y is the Heaviside function. With ^{Et 
—5{Et — E) we found that 



E) = 



Et 



fr 



(51) 



Bo 



which is similar to formula given by dWivanto et a lJ IlQSfjl 
but for an unknown reason with a different numerical factor 
( 3r/(327r) in their case). 

To conclude we could note that the evaporation process 
can also be taken into account by the gaseous model. Follow- 
ing |He£gi^^L|^2^ "^6 can for instance use an additional 
equation based on the velocity escape at radius r law: 

dp r 

-— OC — 

dt ^Et - <i>{r) 

The flux conservation requires that the density distribu- 
tion of t he escape electrons should follow p^sc oc 
iSpitzer a nd Shaoiro 1972). 

For a non isolated case the particles acquire energy 
marginally in excess of escape energy Et — ^{(x) — 
G' Mnie/rt. The energy loss rate is therefore dH = 
—G'niedM/rt and can be seen as the results of the work 
done on the plasma by the external electric field. 



6.4 Stationary solutions 

There is many ways to derive the quasi-stationary solution 
of the Fokker-Planck equation where the left hand side of 
the orbit average Fokker-Planck equation ve rifies df /dt ~ 
so n = Ho is constant. Following iKindFlOGSl we write a first 
order linear differential equation form 

Et 



^ = / 



dE' dE 



f'r'dE' + T 1 f'dE' 
Eq J e 



(52) 



Ho is proportional to the evaporation rate, and A''^ is pro- 
portional to A*' through equation 1501 1. If Do = it is easy 
to fin d that the only solution is a maxwellian one (|Henonl 
Il96(t) which is unphysical because non-truncated. 



6.5 Kramers-Michie-King solutions 

With the same assumption used to derive equation 11611 
namely a thermal bath Maxwellian assumption form for /', 
we could go one step further and simplify the Fokker-Planck 
equation. In the high energy limit {rrieE — £ ^ ksTe) it re- 
duces to an equation, first derived by Kramer (1940) for the 
brownian motion (,Mernikov..l99L) : 



d/ _ rwe(r) d 
dt ^ 2v d£ 



f I df 



(53) 



Our quasi-equilibrium case, ^ ~ 0, immediately leads to 
the Kramer (1940)-Mic hie (1963)-King (1965) type quasi- 
equilibrium jKinJl96d) for the phase-space density function 
/ distribution : 

£ £t 
/oc (e "B-i-i^ -e "BTi^) 
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Figure 4. Comparison between a Monte Carlo simulation, 
the King distribution and the maxwellian one. The left panel 
shows a comparison between the velocity distribution Ne{v) oc 
v''^f(v'^ /2 + ^{ri))dv' taken at the 1/4 shell such as 
NeAfj) = Nel'i. The Mont e Carlo calculation comes from 
iRobicheaux and Hansoi] l2nn.'j figure 1 after 160 ns expansion 
time. The Kramer-King distribution is = 5,r^^ ^ 200 K. The 
right panel shows the radial dependence of the temperature TT 
calculated by formula I56i compare to Ilobicheaux and Hansqn| 
!2003 figure 2 for an expansion time of 125 ns and 500 ns. The 
Maxwellian distribution is radially independent and is clearly not 
adapted. 



We define what we shall call the Kramers-King's tempera- 
ture r<f . 

The electron density ne{r,t) — J^°° f{r,v,t)'inv^dv 
can then be analytically calculated from the King's distri- 
bution: 



*Erf(y^ 




(54) 



where the proportionality factor is obviously n^l/F^lri) 
where is the electron density at the cloud center and 
rj = r^tir — 0) (notation chosen to be identical to the one 
used in BEG evaporations theories). rit{r,t) — ^'^''T^'^jt-t^r**'' 
is solution of the self-consistent Poisson equation: 



r dr 



^{rr]t{r,t)) 



kBT^it) 



qini{r,t) -\- qe,n^(r,t) 



eo 



(55) 



Similarly, with zero velocity of mass transport {kBTe 
meO-v)) we have the important relation 



1 - 



5/2 



ksTr it)Eri{0.22r^t){56) 



with NeTe = 4:TTr'^ne{r)Te(r). We could noticed that, for 
equations 15411 and 1561 1. the adding factor to e^'Erf{^/x) is 
always the power series expansion for e^Erf(-ya;) about the 
point 0. 

At first approximation we can say that electrons are 
always in King's type quasi-equilibrium during the whole 
expansion of the ionic and electronic cloud. Indeed, in fig- 
ure m using Figs. 1 and 2 M onte Carlo simulations of 
IRobicheaux and HansonI (l2003l) . we demonstrate that the 
King electron distribution is better than a pure Maxwellian 
one. 

The harmonic potential approximation in equation 1551 
leads to: 



7 T^-ff ^ ?e 2/0 On 

3eo 



(57) 
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Figure 5. Comparison between numerical calculation and for- 
mula 1581 the square data comes form several A^e , Ni — Ne couples 
and different rj values (3 < ?? < 15) with rt = 15 — 20. 



This equation is in fact well confirmed by numerical reso- 
lution of equation ^ for 10 K< < 1000 K. This ex- 
pression has to be compared with equations JSJ and 12111 
because the electron are not in a gaussian distribution so 
71° / Ne/(2Tva^)'^^^. In fact, we found that ni{r) = ne{r) 
for r/a « 3 — 4 almost independently of all the other pa- 
rameters and that the electron density distribution varies as 
r~^^^ near infinity as expected. These results can be used to 
improved the naive formula extract form the Virial theorem. 
The numerical results can be approximated, see figure |S] by; 



1.9(?7 - 2)fcsr/ 



(58) 
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This is close to the intuitive results: the trapping depth 
•qkuT^^ has to be roughly equal to the trap depth, calcu- 
lated assuming a gaussian shape for ions and electrons where 
1.9(77 — 2) replace the naively rj expected value. 

6.6 Second order solutions 

A simple extension of our work would be to look at the 
evolution of the key parameter 77(t). An equation for the 
sole T) par ameter can be de rived from the entropic variational 
principle jTakahashilll993h : 



: 



Et 



dt 



9 In / dr 



n 



d / dlnf \ 
[ dE ) 



dr) dE dr] 



dE 



(59) 



Another possible equation will be the one based on cor- 
rect mass, e nergy and entropy e volution as done in clus- 
ter physics_jL^inai^m^zgi[ lO^ jfl or in tiltra-c old atomics 
physics l|Berg-SorenseJl997l:lPinke et alJll998l) . 



In the article we have only developed the isotropic 
case. Obviously another simple extension of our work will 
be to use anisotropic model. Similar Fokker-Planck equa- 
tions holds for more general anisotropic system. As a sim- 
ple example we could indicate a natural form for the elec- 
trons quasi-equilibrium distribution: the Michie's one used 
in anisotropic sample: 



f<xe' 



where the second integral J is the angular momentum. It is 
possible to average f{E, J, t) over J to restore the equation 



for fCg) iCohnlll979l) . In such Michie distribution the elec- 



>,3/2 -7/2 



near rt 



tron density distribut i on va ries as {rt 
iBinnev and Tremainel 

To study the general properties o f the plasma it would 
be profitable to go one step further JPratalll97ll) . Indeed 
to study the velocity, or energy, distribution the first order 
is sufficient because these quantities are mainly determined 
by the "maxwellian" center where the collisions occurs. On 
the contrary, the escape rate is also determined by the ve- 
locity distributions in the whole system where / is better 
approximated by the King distribution. 

Using equation 1441 we solve the first order linear dif- 
ferential equation 15211 : 



fiE) 



Ho 



ATvTH^iE') 



r 



'dE' 



(60) 



This gives the new / function and then the news p, "I>, r func- 
tions from equations 1121 and 14511 . Using formulas 1141 and 
lUSl and assuming the Kramers- King's distribution leads to 
10% accuracy (for 1 < 77 and for all possible E values) ap- 
proximation 



TV, 



-'7/2 



{E - Eo) 



(61) 
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The {E - Eo f behavior near the core region E ~ Eo leads 
to a diverging integral in equation 1601 . 

The model looks pathological. But, with the King dis- 
tribution, we can see th at the energy flux is not zero at the 
center. By analogy with lHenonlll96ll we may say that the 
energy flux is due to the Three-Body encounters. The (two 
body) Fokker-Planck equation contain the fact the three 
body recombination does exist ! The energy flux emerging 
from the center is in fact realistic because it is supplied 
by binary formation. Therefore, we need to put the three 
body collisions inside the Fokk er-Planck equation . By anal- 
ogy with the cluster case (e.d. iTakahashU |l99lf)) we must 
add an energy source term to the first-order diffusion coef- 
ficient ijLvman SDitzer^ll987^ (2.79). Using the orbit averag- 
ing of the heating rate HlOjl we need to replace the flux 11 in 
equation ^ by fl = 11 - iV/ with : 



iV = 16n 



-l{i5) 



E{2{E-<i>{r,t)f/\ 



(62) 



The only change is to change A'' by A'' -f Nt, so if N is 
proportional to Nt we recover, through equations H60|l and 
1611 . the King solution. This is a strong indication of the 
good validity of this kind of solution as already proved by 
our comparison with Monte Carlo simulations. In the cluster 
case the number of binary is always negligible and is usually 
not taken into account in this modified Fokker-Planck type 
equation. Equation l|10|l indicates that this assumption is 
not as good in a plasma sample, typically twenty (100/5.4) 
times worse. We then should add a loss term on the right side 
of equation 1481 . A better solution m ight be to use the or- 
bit average m aster equation (s ec Mans bach and Keck! (|W6^ 
(IV. 1) or in iGoodman and Hut. 1,199,'j) (2.21)). The rate 
constant is given by equation with (kBTe)^ K{E, E') = 
mik{m,{E - EQ)/{kBT,),me{E' - £o)/(fcflTe)) : 
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+ 00 

[K{E, E')f(E\ t)^K{E', E)f{E, t)]AE' 

00 

We shall not develop here further on this discus- 
sion but we could mention that the Fokker-Planck equa- 
tions with three body heating terms is similar to the self 
similar Fokker-Planck equa tion without binary heating by 
iHeggie and Stevenson! Il98^ especially with N oc qf and E 
constant. 




7 EVAPORATION EXPERIMENTS AND 
ELECTRON TEMPERATURE 
MEASUREMENT 

One of the most interesting part in the ultra-cold plasma 
system is its ability to be experimentally tested. We have 
already mentioned some experiments as the plasma expan- 
sion (through electron density measurement) or the TBR 
processes. 

Here, we would like to test the phase state distribu- 
tion /. It could be experimentally te sted for instance using 
an electric pulse to extract electrons iVanhaecke et al]l2004l : 
iRoberts et al]l2004) si milarl y to "runaway electron" exper- 
iments (iKvdgrud et alJll973l) . A similar technique has been 
used in neutral atom Bose Einstein C ondensation ( EEC ) 
confined in a static external potential llDovle elTalJilOsjl) . 
The extraction should lead to an instantaneous picture of 
the energy distribution of the electrons in the plasma. But 
our case is complex due to the fact the potential depends on 
the number of trapped particles through Poisson's equation. 

7.1 Experimental determination of the 
temperature. 

Electronic temperature is a key parameter in plasma. If the 
thermal energy is lower than the Coulomb interaction energy 
the plasma approaches the strongly coupled regime where 
correlation effects become important (Pohlctal. 2004c). 
The basic idea we have used in jVanhaecke et alJl2004) was 
to experimentally determine the temperature using a short 
voltage pulse V to extract electrons. The number of elec- 
trons ejected by the voltage V is plotted in Fig. |S|both for 
the plasma only, and for the plasma plus Rydberg sample. 

The first step of the theory is to use the threshold value 
Vi**" (see Fig. [SJ that is necessary to remove all the free elec- 
trons to find n1. Assuming there are always a few electrons 
with zero velocity in the Lagrange point where electrons are 
extracted from the ion potential well, the maximum electric 
field F^^ created by the ionic space charge is: 

F'^-^Y^^ 2M-^n°V2^ (64) 
d 47reo 

Analysis of Fig.|S|and use of formula 1641 lead to knowledge 
of nl at few percent accuracy (assuming a is known exactly). 
Experimentally, one main uncertainty is the imprecise deter- 
mination of A^e due to laser fluctuations and poor calibra- 
tion of the charged particle detector (Micro Channels Plate 
MCP) detector. Our data were interpreted using a Kramer- 
Michie-King electron distribution. The last step of the the- 
ory is based on the solution of Poisson's equation with n° 
known from equation 1641 . Because some electrons are re- 
moved even for small Vi values, rie (r) this indicates that rt ~ 



Figure 6. Number of electrons ejected when varying the voltage 
V applied after 1 fis plasma expansion time. The dashed lines are 
guide for the eye to interpret equation I64i . Each data have been 
calibrated to the average value of the total electron number. 



ex. In order to take into account our magnetic and stray elec- 
tric field we shall assume rt « 12a. Furthermore, the total 
number of calculated electrons Ne{t) — J^°° ne{r,t)i'Kr^dr 
should reproduce the observed number A^'e. This determines 
the values for the two remaining parameters n° and as a 
function of the unknown rj parameter through formula 15811 . 
Unfortunately, we were too sensitive to the exact value of 
to be able to give an absolute v alue for the plasma temper - 
ature. Figure |I| and Fig. 10 of jKuzmin and O'Neill20023) 
show that rj 8 seems a reasonable choice to interpret the 
data given in figure |S| Even with this not fully sat isfactory 
assumption we have used in IVanhaecke et alj(2004l the the- 
ory to give relative results. One conclusion was that the 
temperature should not increase or decrease by more than 
a factor 5 when Rydberg atoms are added into a plasma. 

The theory presented here is based on several assump- 
tion which have to be tested. First, to check the validity 
of formula 164t we have compared, after 2 /is of expansion 
time, the plasma size obtained by formula 1641 to the one 
given by the expression ^ for different plasma expansion 
velocity vq. Using = ^2-^3/2 1 equation 1641 leads to 

a oc NiVi^. Ni is known from MCP The experimental re- 
sults are depicted in figure|7|where Vq ~ ksT^ /rrii is varied 
by changing the laser ionization wavelength. Fitting results 
are in agreement with formula JTJ for ctq ~ 200 /im which 
is very close to the expected MOT size. Therefore, we have 
here a strong indication concerning the validity of formula 

Second, we have compared formula 1581 . which gives 
, with . has been experimentally extracted using 
formula 1581 . with 77 = 7, cr is obtained using the results in 
figure |7| Ni is given by equation 1581 knowing the voltage 
threshold value Vl^ . Results in figure JHJ, for three different 
ion numbers, seems to indicate that is close to and 
that = 7 is a reasonable value, after 2 of expansion 
time, but has probably to be adapted for each ion number 
case. 

Up to now we have just used one single point of the 
full curve given in figure |S] a more complete study should 
yield to a determination of the r] value. The full process is 
beyond the scope of the article but might be tackled with 
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Figure 7. Plasma size extract from expression ^ = ■^^ir 

compared to ^ k^T2 jriXi wliicli is known from the laser ion- 
ization wavelength. a^^N^ are respectively the sample size 
(assumed to be frozen), the ion number and the threshold voltage 
for the laser sets at the ionization threshold (Tj « 0). cr,Ni, V^^ 
are respectively the sample size, the ion number and the threshold 
voltage for different laser wavelength taken after 2 /xs of expansion 
time. The linear fit is cr^ CTq + 10~*Dg. 
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Figure 8. Kramers-King temperature , extracts from equa- 
tion HHJ (with = 7), versus laser wavelength, (a) The photoion- 
ization electronic temperature T2 defined by: 'ikBT2 /2 equals 
to the difference between the photon energy and the ionization 
threshold (near 508 nm). (b) data taken with full laser power 
(large A''^ value), (c) with 5 times lower power than in (b) (medium 
Ni number, and (d) with 15 times lower power than in (b) (small 
Ni value). 



the "violent relaxation" theory llZiegle r and WiecherJll989t 
IWiechen and Zieglejll994l: IChavanisl ll998). Indeed, even m 
an impulse approximation with Fokker-Planck-Vlasov col- 
lisionless equations, or adiabatic approximation in gaseous 
equations, the violent ejection is a very complex system. At 
the present experimental status we have only investigated a 
global parameter namely the electronic temperature of the 
plasma. 



7.2 Other methods 

Many other methods may be uses to determined the plasma 
temperature. 

• IWalz-Flannigan et al.ll2004l have compared an experi- 
mental Rydberg binding energy distribution with calculated 
outcomes of inelastic collisions between Rydberg atoms and 
electrons in the plasma taken from equations @. They con- 
clude that it is very likely that electron-Rydberg atom col- 
lisions cause most of the experimentally observed popula- 
tion redistribution into states. Their fitting method can then 
serve as a tool to estimate the electron temperature. 

• The plasma expansion itself is already a signature of the 
(time average) temperature as indicated by equations 131)1 
and l|7|l. One pr omising way toward this measurement has 
been provided bv lKillian et al .*2003'who use Strontium ions. 
Indeed, the non alkali ions as the strontium one have visi- 
ble resonance transitions and can then be detected by laser 
induced fluorescence. The Doppler shift of the resonance in- 
dicates that the ion velocity, and so it's vari ation (the ac- 
celeration), reflects the electron temperature iSimien et alJ 
12004) . 

For alkali atoms it is no more possible to use optical de- 
tection. We might think to use position sensitive detector to 
image the ionic cloud. A similar idea is to use the time of 
flight projection image of the ionic cloud. We have tried this 
technique. It is efflcient and easy to handle only for really 
small ion number clouds to avoid a too strong Coulomb ex- 
plosion during the ion transport toward the MCP. Indeed, 
even using formula I32L which is based on a hypothetical 
pure spherical symmetry, it is hard to restore the initial ionic 
cloud shape from the projection image of the ions in the de- 
tector. 

• Another method to study the plasma temperature is to 
monitor the electron evaporation. Experimentally, the eas- 
iest way is to add a small electric field to force slow evap- 
oration of particles. Theoretically this is the well known 
Kramers problem of escape from a trap (lMernikovlll99ll: 
iPeter et alJll99(] h. The escape rate (see equation (I51II ). can 
be calculated for instance using a Ki ng distribution and is 
usual ly found to be close to l/(100io) iBinnev and Tremaind 
Il987l) where te = 9al / [IQ-JlxG''^ mlnlhi K) is the electron 
thermalisation time. This results assume a purely energetic 
truncated condition. But, in globular cluster dynamics it is 
well known that this hypothesis of purely energetic trun- 
cated con dition is not wel l suite d especially for aniso tropic 
systems jTakahashi et alj Il997l: iKim and Obi Il999t) . The 
comparison between ultra-cold plasma experiments and the 
theories might help to understand how accurate the assump- 
tion of energy truncation is. This technique has been pro- 
posed by the NIST group using a constant small electric 
field to detect all the eject electrons. Some results are re- 
ported in Robichcaux and Hanson 20q3 figure 3 but were 
not interpreted. The results are in fact compatible with the 
l/(100fc) evaporation rate on the order of lO'^^ s~^ decreas- 
ing with time due to the fact the density decrease much 
faster than the electron average velocity during the plasma 
expansion. 

• Similarly to our work, iRoberts et aljl2'004l used electric 
field pulse to extract the plasma temperature. To avoid any 
discussion concerning the potential shape they tried to affect 
electrons only in the asymptotic part (in 1/r) of the trap- 
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ping potential bu using only very small electric fields. Their 
nice experimental data are treated using a theory based on 
maxwellian electron distribution which we believed is not 
as appropriate as the Kramers-King's one, especially for the 
slightly bound electrons which have been removed by the 
electric pulse. Furthermore, even with a small electron frac- 
tion removed, A'^; — TVe can be affected and the electron ex- 
traction process is probably very complex. The main results 
concerns the electron cooling during the plasma expansion 
and the almost constant initial temperature ~ 50K in- 
dependently of Tg . This final results looks surprising to us 
because it seems to disagree with the Virial based equations 
12111 15711 or 1581 1. Unfortunately, the published data are not 
sufficient to check the results with other electron distribu- 
tion. 



8 CONCLUSION 

We have shown that the dynamical behavior of the electrons 
in an ultracold plasma is similar to the one described by con- 
ventional models of stars clusters dynamics. The evolution 
of such a sample is dominated by kinetics equations formally 
identical to the ones controlling the evolution of a stars clus- 
ter using only a new negative gravitational constant defined 
by equation Q. We have developed here some aspects of the 
analogy with globular star cluster dynamics. The infiuence 
of stellar encounters namely: relaxation, equipartition, mass 
segregation, escape, inelastic encounters (coalescence, dissi- 
pation of energy), binary formation by three body encoun- 
ters, interaction with primordial binaries stars,... have exact 
partners in the plasma case. There is similar law for the three 
body recombination for instance the Heggie's law: "hard bi- 
naries get harder and soft binaries get softer" is identical in 
both system. Thanks to the Virial theorem we could relate 
the plasma temperature with the ions and electrons num- 
bers (or central density) and with the initial kinetic energy 
given to electrons by the laser ionization. We found that the 
Fokker-Planck equation is exactly identical for electron in 
ultra cold plasma and for stars in globular cluster. The only 
modification occurs in the potential which, in the plasma 
case, should contain an extra part due to the ions. The ions 
are spectator but are needed to create a confining potential 
for the electrons. The gaseous equations, coming from the 
first moments equations of the Fokker-Planck equation, con- 
firm the early self similar evolution of the ionic cloud. This 
self similar evolution breakdown after few microsecond and 
we have analytically studied the ion shock wave occurring 
when the quasineutrality is violated. The electrons density 
evolution can be studied using similar ffuid equations than 
the one used for clusters. The orbit average study of the 
Fokker-Planck equation indicates that the quasi-static solu- 
tion for the phase space density electron distribution / is a 
Kramers-King's one / oc e^^^^^^' ' — 1 where deals for 
King (or Kramers) temperature although strictly speaking 
a thermodynamics temperature is not defined for a nonequi- 
librium distribution. This is confirmed by comparison with 
Monte Carlo simulation previously published. Even when 
the three body collisions are taken into account this approx- 
imate distribution seems accurate. In fact several numerical 
codes developed for stellar system could be easily adapted 
to treat the ultracold plasma case. Finally, it is possible to 



experimentally simulate the tidal escape from globular clus- 
ter under external galaxy gravitational attraction by using 
small external field to extract electron trapped in an ionic 
cloud. The full process is beyond the scope of the article 
but might be tackled with the violent relaxation theory. We 
have only study here the link between the electrons and ions 
numbers, the plasma temperature and the charged particles 
potential depth through formula 1581 1. 

We hope that this article will stimulate links with the 
astrophysics community and we hope the analogy is also 
useful to develop more physical insight on the ultra-cold 
plas ma physics behavio r. For instance as in the cluster 
case iMevl an and Heggid 1^^) binary collisions might be 
a route for the formation of long-lived multiple systems 
as three-body stable system s or for the forma tion of gi- 
ant molecules as "trilobites" (iGreene et al or macro- 
Rydberg ones in the plasma case iFaroog^^Li SOOS'l. The 
ultracold plasma experiments evolve rapidly, we have men- 
tione d the laser m anipulation of the ions in the strontium 
case iSimien et al ]|2004). This might lead to control of the 
ion motions and to cool sample where crystallization might 
appear (Pohl et ah 2004c). This manipulation of ions can be 
added with external field manipulation as RF electric field, 
magnetic gradient manipulation and might be able to sim- 
ulate a lot of events as tidal shock or tidal heating in clus- 
ter physics. Some picture s of the ionic cloud have already 
been done bv lSimien et al. 12004) and this gives exactly the 
projected density profile in strong analogy with the main 
observed quantity in star cluster: the projected mass pro- 
file. A charge particle detector (MCP) with position sensi- 
tivity might also give similar data. Therefore, analogy with 
dynamics of globular clusters might stimulate new ideas in 
ultra-cold plasma physics. One advantage of plasma system 
compared to globular cluster one results in its physical be- 
havior very closely related to what can be simulated. But, 
the major advantage is the experimental capability of tuning 
different parameters in a huge range. In globular cluster, non 
evolutionary models or single mass models are often used 
but are far from reality, this is no more the case in ultra- 
cold plasma physics. We believe that a comparison between 
different models with very well controlled experiments will 
help to improve further development. 
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